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THE DISTRIBUTION OF BAROCLINITY UITHIN THE ATMOSPHERE 


D.A. Barber and M. Wai, Dept, of Atmospheric Sciences, Oregon 
State University. 


ABSTRACT 


A three dimensional numerical interpolation scheme which 
resolves frontal gradients with high fidelity has been 
developed and is being applied to the study of atmospheric 
upper baroclinic zones. 

INTRODUCTION. 

The aim of our research is to impro/e understanding of the 
structure and evolution of atmospheric fronts. Achievement of 
this goal requires the representation of frontal intensity with a 
dynamically significant parameter for which we have selected the 
baroclinity. 
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where ^ is the baroclinity vector, g the acceleration of gravity, 
and 6 and Vp0 are the potential temperature and its isobaric gra- 
dient. Furthermore, we may define local frontogenesis to be the 
local tendency of the baroclinity which may be shown to be (to 
good approximation) 
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Because of the small horizontal scale of fronts compared with 
the spacing of sounding stations the magnitude of the baroclinity 
may be considerably underest imatea in the vicinity of intense 
fronts if only horizontal or isobaric data are utilized. This 
difficulty may be largely obviated through use of cross-section 
analyses which effectively apply the high resolution information 
available in the vertical over eacn sounding station to improve- 
ment of horizontal resolution. (Shapiro, 1970). I 
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The considerable expenditure of effort required for manual 
analysis of cross-sections has led to the development of success- 
ful numerical interpolation schemes for cross-section analysis 
(Shapiro and Hastings, 1973 and Whittaker and Petersen, 1975); 
however, these are two dimensional techniques from which one can 
obtain isobaric analyses only by laborious combination of many 
independently analysed cross-sections. We propose instead to 
apply a fully three dimensional scheme based on a straight forward 
extension of the Barnes (1973) horizontal method. 


We extend the Barnes scheme by redefining his weight as 



where W^ is the weight factor for the i^^ observation, d^ and Zj 

are the horizontal and vertical distances between the observation 
and a given grid point, and k and D are weight parameters con- 
trolling the smoothness of the output. To facilitate selection of 
k and D, we have coded a two dimensional cross-section analysis 
method based on (3) having a horizontal grid spacing of one lati- 
tude degree (111 km) and a vertical spacing of 50 mb. 

The resulting potential temperature analysis on a cross 
section through the complex hyperbaroclinic zone studied by Frank 
and Barber (1977) is shown in Fig. 1 together with the component 
of the baroclinity in the plane of the cross-section estimated by 
finite difference method. For comparison, Fig. 2 reproduces the 
careful manual analysis of Frank and Barber. The numerical scheme 
has reproduced not only the gross features, but also most of the 
details with considerable fidelity. Further improvement would be 
possible (at the expense of more computer time) if the vertical 
resolution were increased. 

A quantitative demonstration of the efficacy of the cross- 
section scheme is illustrated in Fig. 3. The baroclinity estimated 
from the cross section at the 500 and 800 mb levels is shown by 
the dashed lines. The solid lines are the results obtained 
through application of a horizontal Barnes analysis on a one 
latitude degree by one longitude degree grid. The magnitude of 
the baroclinity estimated from the cross-section program is fifty 
to one hundred percent greater in the hyperbaroclinic zones than 
that obtained from the isobaric analyses. 

CURRENT WORK ; 

We are now coding the full three dimensional scheme for appli- 
cation to the AVF. II data. Furthermore, we intend to demonstrate : 
its usefulness through application to operationally available i 



Fig. 1. Objectively analyzed cross-section through hyper- 
baroclinic zone. Isentropes (<) solid and baroclinity 
(10-'^ s‘2) dashed. 1115 uT, 12 May 1974. 
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Fig. 2. Manually analyzed cross-section along same line as 
Fig. 1. Isentropes (<) solid and isotachs (<t) dashed of 
flow normal to section (Frank and Barber, 1977). 
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Fig. 3. Baroclinity s"^) from objective isobaric analysis 

(solid) and from objective cross-section (cashed) as a function 
of latitude along line of Fig. 1. 
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teletype data and comparison with National Weather Service analyses 

In addition to depicting the full three dimensional baro- 
clinity distribution, wo plan to apply the method to the wind 
field. In the case of the wind field, some smoothing of the 
vertical wind field is desirable (see Schmidt and Johnson, 1972). 
Therefore, a different set of weight parameters (k and D) will be 
appropriate. 

Finally, we shall apply the improved three dimensional wind 
and temperature analyses to the evaluation of terms in the 
frontogenesis equation (2) . 
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RECENT PROGRESS ON WORK UNDER CONTRACT NAS 8-32694 HAS INCLUDED THE FOLLOWING 
MAJOR ITEMS. 

% 

1. A preliminary teit of the full three dimensional numerical interpolation 
scheme has been applied to the wind field for OOOO UT, 12 May 1974, over a 
limited region of the midwest. The area is restricted so that testing of the 
scheme may be accomplished without undue expenditure of computer time. The 
results are shown in Figures 1-7 a, b and c. These figures show the wind speed 
and direction (a), the isobaric divergence (b) , and the isobaric vorticity (c) 

ft 

for 7 selected pressure levels on a 1 latitud>' degree by 1 longitude degree 
grid. 

With two important exceptions, these figures demonstrate excellent vertical 
and horizontal continuity and quite reasonable fields of the kinematic variables. 
We believe that, when applied to a larger sample of data, the method will prove 
a significant improvement over existing techniques. 

The two exceptions noted above are related to horizontal and vertical 
boundary |>."'*'l'jms. The first of these will be substantially overcome by merely 
including a larger area within the analysis domain, ".'he second results from 
difficulty near the surface of the earth. This problem was encountered in 
earlier tests of the scheme for analysis of the potential temperature field; 
it was overcome by including bogus data points below the earth's surface obtained 
by downward extrapolation based on data above the surfice. This device has not 
(icon successful in the wind case. 

2 . We have obtained a small grant of computer time from the Oregon State 
University Computer Center for the testing of the three dimensional scheme on 
u case involving use of standard teletype transmissions of sounding data. The 
temperature data for 00 UT, 2(> January 1978, have been punched for this purpose. 
While this work is not directly supported by NAS 8-32694 it will demonstrate 


- 2 - 


the usefulness of the three dimensional scheme applied to conventional, 
operationally available data. 

3. The often frustrating efforts to convert the data tape for the every 

contact AVC II rawinsonde data to a format comba table with our local COC 

computer have continued. The OSU Computer Center has been unable to perform 

this task so we have sent the tape elsewhere for conversion. 

Immediate future work will emphasize the following items: 

1. Improvcmcnc of the bottom boundary condition for the analysis scheme. 

2. Preparation of programs to compute vertical velocity by the kinematic 

method iicluding standard correction procedures. 

3. Enlarging the area for testing the analysis scheme and adding 

temperature data to the wind data for the test case. 

4. Testing schemes to reduce the running ..irac of the analysis scheme. 

5. Preparation to run analyses on several time periods of AVE II data 

when the data tape has been converted. 

6. Running the analysis program on the 26 January 1978 case to illustrate 

the usefulness of the technique when applied to conventional data. 
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